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We discuss the correlators in lattice HQET that are needed to go beyond the static theory. Based on our 
implementation in the Schrodinger functional we focus on their signal-to-noise ratios and check that a reasonable 
statistical precision can be reached in quantities like /b 3 and Mb* — Mb- 



1. INTRODUCTION 

The physics of mixings and decays of B-mesons 
is crucial for a determination of several entries of 
the CKM-matrix, e.g. V u b- To relate experimen- 
tal observations to Standard Model parameters, 
transition elements of the effective weak Hamil- 
tonian must be computed in a reliable fully non- 
perturbative framework, e.g. on the lattice. 

Since the large mass of the 6-quark evades a 
direct treatment at lattice spacings a -1 = 2 — 
4 GeV, effective field theory methods are invoked. 
The one tried first is the static approximation 
which is the leading order of an expansion in 
the inverse heavy mass, i.e. the leading order of 
HQET £Q. Here, the heavy flavor field ifa satisfies 
(1 + 7o)V'ii = 2V>h and has a discretized action 

sr^^^W^hW (i) 

X 

D iph{x) = [iph(x)-U^(x-a6,0)i> b (x-a6)]/a . 

A non-perturbative renormalization of the effec- 
tive theory - including the subleading terms - 
is possible, i.e. divergences can be subtracted im- 
plicitly through a direct matching against QCD in 
a small volume, and very precise results at short 
(xq < 0.5 fm) distances may be obtained 2 . One 
of the advantages of such a strategy is that the 
theory maintains a well-defined continuum limit, 
while this would not be true, if the matching was 
performed to any fixed order in perturbation the- 
ory. We consider this important. 



There are, however, some intrinsic difficulties 
with this approach. The first one is the statis- 
tical noise in correlators at large (xq > 0.5 fm) 
euclidean distances. Brute force is no solution, 
since the noise- over- signal ratio grows exponen- 
tially in xq. This is related to the 1/a divergence 
that has been subtracted. If the HQET approach 
is extended beyond the leading order, l/a™ diver- 
gences are implicitly canceled, and the obvious 
fear is that the noise problem will be even worse. 
The second difficulty relates to the fact that any 
limited precision in the matching will translate 
into additional uncertainties of all quantities in 
the large volume; hence devising good matching 
conditions is an important task. Here, we shall 
address point one, i.e. the statistical precision 
that can be reached in typical correlators at the 
next-to-leading order in the HQET expansion. 

2. CORRELATION FUNCTIONS 

The Schrodinger functional (SF) master for- 
mula for the heavy meson decay constant 1 
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relates the expression in the relativistic formula- 
tion to that in the HQET approach. In the latter, 
the correlation functions at leading order are 
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A stat = -^(Q 'o) 

1 For any unexplained notation we refer to 12131 . 



(2) 
(3) 



1 



2 



with the static-light axial current 
Af at {x) = $(x)~/ o75 ip h (x) 
and the boundary operators 

O = VCh(y)75^(y-z)C(z) 



(4) 



O' = 



E 

y',z' 



C'(y')75^(y'-z')^K) 



on the bottom and top (x — 0, T) of the SF-box. 

The idea of lattice HQET at subleading order is 
to expand the Lagrangian (we include coefficients 
for subsequent renormalization) 
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2m 2m 
in the exponent and to treat the new terms as 
insertions in any correlator. Expanding consis- 
tently means that one keeps just terms in 1/m, 
no products 0(l/m 2 ). The correlator f A ta,t is thus 
augmented by a) kln / spin /(2m) times 

l « at (x) V A kin/spi »0) (5) 
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with X km (u) = ■0h(u)D 2 '(/'h(u) as bulk insertion 
or X s P in (u)=if> h (u)a-BiP h (u), and ditto for /f at . 
The dimension 5 pieces of the Lagrangian appear 
only as insertions, and this is crucial for the renor- 
malizability [to order 1 jm\ of the theory. 
The expansion of /a to order 1/m reads 
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where w kln = u) kln /(2m) (= ^oa, in the notation 
of [2]) and similarly for w spin . An analogous ex- 
pression (without the <5/ A tat piece defined in|3]) 
replaces /f at . The coefficients a (o/i) ;W kin/s P i n 
are functions of M;,o, g$ and may be determined 
as in [2] • Hence 
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+ (6) 
u> spin R spin (T/2)} 



where 
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In the same way, an effective mass m e s(xo) - 
a^ 1 log(fji L (xo)/ fA(xg+a)) is 1/m-expanded as 

TOeff(zo) = m^ t (x )+Sm + 



LO'""r 

with mfg t (xo) from / A 



Jdn (x ) +uj spiD r spin (x ) (7) 
stat (a;o) and 



r kln (x ) = 
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Here, the contribution 



/ A tat (^ + Q 



has been suppressed, since in the plateau region 
the two terms would cancel, while r km / s P m ap- 
proach a constant value each. 

3. NUMERICAL RESULTS 

The problem of large fluctuations in the static- 
light propagator was solved by modifying the co- 
variant derivative in QJ. As shown in 3 , an 
APE- or HYP-smeared link V^{x — aO, 0) instead 
of U' leads to an improvement in the noise-over- 
signal ratio of / A tat (:co) which grows exponen- 
tially with Xq while the discretization errors re- 
main at the same level. Fig. ^ shows that a 
similar improvement is found at subleading or- 
der in the 1/m expansion. HYP-links 4 with 
(ai, CK2, CX3) = (1.0,1.0,0.5) do better than those 
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Figure 1. Noise-over-signal ratio for the heavy- 
light correlator / kln , using three discretizations 
of the heavy propagator. Similar graphs are ob- 
tained at leading order in the 1 jm expansion 3 . 
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with (ai, a-z, a^) = (0.75, 0.6, 0.3), which, in turn, 
are better than APE-links 5 with staple-only 
contribution and no S'?7(3)-projection. The same 
ordering was found at leading order. 

At this point we cannot compute the sublead- 
ing contribution to /b b , since the matching co- 
efficients a^^ x \ cj km / s P m have not yet been de- 
termined, but we can assess the statistical error 
of such a contribution. Our figure and numbers 
stem from a quenched simulation in a 16 3 x 24 
SF-box at (3 — 6.0 with k ~ ^strange and 5600 
measurements, and we shall quote only the re- 
sults for the improved HYP-parametrization. For 
the two correction terms in JSJ we find the values 
ai? kin = -0.78(14), aR B P in = -0.83(03). 

For an interpretation of their statistical er- 
rors, we need a rough figure for the coeffi- 
cients w km / s P m . Up to typical (logarithmically in- 
dependent) renormalization factors of order one, 
they can be estimated by their tree-level values 
a -i w ki„/s P in _ i/(2 a m) ~ 0.2 (from a" 1 ~ 2 GeV 

and m~5 GeV). In addition, the renormalization 
of i? kin /spin does 

require power divergent subtrac- 
tions, but they are contained in (for which, 
therefore, a tree-level estimate would be no good). 

From kj km / s P m ~ 0.2 one gets a statistical er- 
ror of the kinetic contribution of 3% and an er- 
ror of the spin contribution of 0.6% in expression 
©■ Neglecting uncertainties in the c^ / 1 ) and 
cj km / spm , it appears that already with the meth- 
ods applied here, one can compute /b s at sub- 
leading order in the HQET expansion in such a 
way that the noise in the dimension 5 correlators 
increases the error by ^3-4 percentage points. 

A simpler application, where even a first es- 
timate may be given, is the vector-pseudoscalar 
splitting Mb*— Mb- This effect sets in at sublead- 
ing order in the HQET expansion, 

Mb* — Mb = ^w spin r spin (a;o) (8) 

and has first been estimated in a lattice computa- 
tion in . Our result for the large-^o asymptotic 
is a 2 r spin = 0.060(5). With a~V pin ~ 0.2 one 
gets M B * -M b = 0.016(1) (2GeV) =0.032(3) GeV, 
which should be compared to the experimental 
result 0.046 GeV. Hence, the lattice value falls 
short by about a third. Whether this reflects a 



large renormalization factor (similar to csw)j a 
slow convergence of the HQET expansion or a 
genuine quenching artefact remains to be seen, al- 
though the results of [7] render the second expla- 
nation rather unlikely. We emphasize that these 
figures are just indicative; reliable numbers can be 
given only when to spin has been determined and 
smaller lattice spacings have been considered. 

4. CONCLUSION 

We have attempted a first test with correla- 
tor ratios needed to compute observables in the 
heavy-light system beyond the static approxima- 
tion. The natural fear that they will, for large eu- 
clidean distances, be even noisier than the leading 
order piece seems not to become true. This is be- 
cause the reduction of the noise-over-signal ratio 
via a better HQ-discretization works more effi- 
ciently at subleading order. The statistical preci- 
sion may be further improved, e.g. by employing 
a smaller time extent in ff at and j kln / spm (to- 
gether with a careful choice of the wave function, 
see [3]) or by the methods of 0. We conclude that 
lattice HQET at subleading order looks promising 
enough to motivate further studies. 

We thank A. Shindler and M. Delia Morte for 
checks on the implementation and discussions. 
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